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Abstract 

An  engineering  approach  to  the  design  of  an  in-line  dynamic  absorber  is  presented  using 
transfer  matrix  techniques  for  an  infinitely  long  cable  system.  The  natural  frequencies  of 
the  system  are  determined  and  the  properties  of  the  installed  damper  are  optimized  to 
maximize  the  amount  of  energy  absorbed.  Numerical  simulations  show  that  the  optimally 
tuned  and  damped  in-line  dynamic  absorber  will  dissipate  50%  of  incident  transverse 
vibration  energy.  A  design  methodology  for  a  nominal  marine  monitoring  program  is 
presented  and  specific  engineering  constraints  provided  to  ensure  that  the  modeled  system 
meets  design  specifications 
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Chapter  1 


Introduction 

Investigation  of  the  marine  environment  is  a  growing  industry.  In  order  to 
economically  obtain  enough  information,  increased  usage  of  remote  sensing  technologies 
will  be  required.  The  use  of  these  technologies  will  require  that  such  information 
collection  systems  be  robust  enough  to  survive  in  the  harsh  environment  for  prolonged 
periods  of  time  without  human  intervention.    In  most  cases,  the  areas  of  concern  for  the 
scientist  can  be  defined  by  a  geographical  boundary,  and  an  ocean  mooring  can  be  used  to 
locate  the  measuring  devices.  The  motions  of  the  moored  body  can  damage  the 
equipment  or  produce  erroneous  readings.    Thus,  the  robust  requirement  of  the 
measurement  system  will  lead  to  either  a  hardening  of  the  sensitive  electronics  of  the 
system  or  other  some  other  sort  of  isolation  method  that  will  decrease  the  cyclic  loading 
that  will  be  seen  by  the  system. 

The  problem  of  the  motion  of  a  moored  buoy  with  sensitive  electronic  instruments 
located  along  the  length  of  the  cable  can  be  modeled  as  a  string  with  a  rigid  lumps.  A 
previous  analytical  approach  of  the  problem  was  performed  by  Li  [1].  In  his  analysis,  Li 
was  able  to  provide  a  analytical  description  of  how  to  model  a  string/rigid  lump  system  by 
the  geometric  and  material  properties  of  the  components  and  the  cable  and  consider  the 
system  more  than  just  a  lumped  mass  at  a  particular  point  along  the  cable.  The  energy  of 
the  wave  traveling  along  the  cable  might  be  transmitted  through  the  component  or 
reflected  back,  depending  upon  the  specifics  of  the  problem. 

If  the  physical  connection  of  the  equipment  along  the  length  of  the  cable  were  to 
include  a  spring  and  a  damper,  the  model  becomes  more  complex.  The  energy  of  the 
wave  traveling  along  the  cable  might  then  be  absorbed,  transmitted  or  reflected  at  the 
interface,  dependent  upon  the  nature  of  the  system  and  the  frequency  of  the  incoming 
wave.  Li  was  able  to  investigate  the  modeling  of  a  mass/string/damper  system  that  would 


be  installed  along  the  line,  as  well  as  endpoint  terminations  of  the  cable — either  connected 
to  a  moving  or  a  fixed  frame  of  reference.  The  in-line  system  has  more  of  a  practical 
application  as  it  can  become  a  part  of  the  device  itself,  as  components  may  be  changed 
along  the  length  of  a  moored  cable  in  order  to  alter  the  measurement  systems  available  to 
the  scientist  in  the  field. 

The  idea  that  the  in-line  component  could  be  isolated  via  an  appropriately 
engineered  system  has  significant  applications.  The  component  then  could  be  designed 
more  towards  the  gathering  of  information,  not  just  being  able  to  withstand  the  normal 
marine  forces.  Using  the  ideas  established  by  Li,  this  thesis  will  investigate  the 
characteristics  of  such  a  moored  body  system  and  develop  an  engineering  approach  to 
determining  the  optimum  value  of  the  constituent  parts  of  the  measurement  devices  such 
that  the  operation  of  the  in-line  components  is  not  impaired. 

This  thesis  will  investigate  the  design  of  an  in-line  dynamic  absorber  tuned  to  the 
frequency  where  the  maximum  amount  of  energy  traveling  along  the  ocean  mooring 
system  can  be  absorbed.  Chapter  2  outlines  the  rationale  that  drives  the  need  for  this 
analysis.  Additionally,  an  application  that  presents  a  process  for  the  design  of  a  marine 
monitoring  program  and  establishes  technical  requirements  is  introduced. 

Chapter  3  is  the  introduction  to  the  use  of  transfer  matrix  notation  and  the  benefits 
that  result  in  the  use  of  this  system.  Chapter  4  is  the  analysis  section  of  this  thesis.  The 
search  for  an  analytical  approach  to  determining  the  optimum  amount  of  energy  that  can 
be  absorbed  in  the  mooring  system  is  provided.  A  numerical  result  is  obtained  and  the 
optimum  value  for  the  absorber  for  a  given  system  is  provided.  Chapter  5  is  a  discussion 
of  the  major  findings  and  potential  areas  for  further  research. 


Chapter  2 

Oceanographic  Monitoring  System 

2.1       Introduction 

The  design  process  to  implement  an  oceanographic  monitoring  system  discussed  in 
chapters  1  must  minimize  the  effects  of  marine  induced  forces  on  an  infinitely  long  cable 
system  with  installed  monitoring  equipment.  One  practical  application  of  the  design 
process  will  involve  the  use  of  a  monitoring  system  that  is  located  in  a  marine  environment 
of  interest.  This  chapter  will  describe  how  such  a  monitoring  plan  is  designed  and 
implemented. 

The  first  step  in  the  development  of  a  monitoring  plan  requires  an  executive  level 
design  methodology.  Using  the  concepts  from  a  report  by  the  National  Research  Council 
[2],  the  following  steps  are  used: 

1 .  Define  Expectations  and  Goals 

2.  Define  Study  Strategy 

3.  Conduct  Exploratory  Studies  if  Needed 

4.  Develop  Sampling  Design 

5.  Implement  Study 

The  expectations  and  goals  of  the  particular  monitoring  plan  would  be  based  upon 
the  marine  site  chosen.  The  ultimate  goal  is  to  have  a  monitoring  system  that  produces 
useful  information  that  can  used  to  make  a  reasoned  decision.  Thus,  the  effects  that  are 
seen  in  the  monitoring  results  must  reflect  information  about  the  situation  that  exists  in  the 
marine  environment.  The  legal  framework  that  guides  the  administration  of  human  actions 
in  the  marine  environment  also  molds  the  expectations  and  goals  of  the  monitoring  plan. 

Defining  the  study  strategy  of  the  monitoring  plan  requires  the  implementation  of  a 
strategy  that  begins  with  the  identification  of  resources  in  the  marine  environment  that  are 
at  risk.  From  that  identification,  a  conceptual  model  is  developed  that  establishes  links 
that  exist  in  the  marine  ecosystem  components.  These  links  involve  spatial,  temporal, 
biological  and  chemical  variables  that  become  the  basis  of  the  study  strategy.  Feedback 


among  the  monitored  variables  in  the  plan  might  dictate  that  some  be  changed  to  better 
reflect  the  marine  environment.  The  inherent  uncertainties  of  a  conceptual  model  must 
also  be  considered  in  the  definition  of  the  study  strategy.  The  final  result  is  a  set  of 
questions  concerning  the  marine  environment  that  can  be  tested  to  describe  the  state  of  the 
marine  environment. 

Exploratory  studies  that  have  been  performed  on  similar  marine  ecosystems  can  be 
used  as  a  basis  for  the  monitoring  plan.  If  the  need  exists  to  validate  certain  technical 
issues  such  as  the  measurement  of  parameters  under  various  hydrodynamic  conditions  (for 
a  coastal  system,  high  river  flow  from  spring  melt  of  snow  or  low  river  flow  from  a  period 
of  no  precipitation  ),  preliminary  studies  should  be  performed.  The  primary  research 
performed  for  this  validation  becomes  an  interdependent  part  of  the  monitoring  plan, 
altered  as  the  need  arises  to  obtain  specific  information. 

The  development  of  a  sampling  plan  uses  the  testable  questions  described 
previously  in  the  study  plan  and  identifies  meaningful  levels  of  changes  that  can  be 
measured.  In  the  establishment  of  these  meaningful  levels  of  change,  the  natural  variability 
of  the  marine  environment  in  question  needs  to  be  considered.  A  meaningful  level  of 
change  must  also  be  established  that  will  be  observable  in  a  time  frame  that  allows  a 
reasoned  decision  to  be  made  that  will  have  an  impact  on  the  system  being  measured.  This 
definition  of  meaningful  change  will  also  be  as  fluid  as  the  improvement  of  prediction 
models  and  the  analysis  of  data  will  allow.   The  selection  of  the  variables  to  measure  in  the 
monitoring  plan  is  limited  by  the  resources  available,  and  they  must  be  chosen  to  provide 
the  most  useful  information  about  the  status  of  the  system.  Thus,  variables  may  be 
measured  that  will  provide  a  realistic  view  of  the  system,  although  they  may  not  be  part  of 
the  predictive  model.  These  substitute  variables  may  contain  early  warning  signs  of 
system  degradation  or  information  concerning  specific  items  of  interest  to  the  users  of  the 
sampling  plan.  In  either  case,  the  substitute  variables  will  provide  information  that  can  be 
applied  to  the  predictive  model  via  knowledge  of  the  interactions  between  the  substitute 
variables  and  the  model  variables.  The  use  of  specific  oceanographic  measurement 
equipment  is  then  chosen  to  reflect  the  variables  measured. 
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Two  concerns  that  must  be  addressed  in  the  design  of  the  sampling  plan  are  the 
level  of  confidence  in  the  knowledge  of  the  interactions  among  the  elements  of  the  system 
and  the  statistical  characteristics  of  the  sampling  of  the  designated  variables.  If  the 
interactions  among  the  biological  elements  is  purely  theory,  then  the  uncertainty  must  be 
addressed  in  the  analysis.  If  the  sampling  process  requires  a  significant  number  of  samples 
in  order  to  verify  the  validity  of  the  measurement,  then  the  monitoring  plan  is  not  efficient 
in  obtaining  the  maximum  amount  of  information  for  the  effort. 

Using  the  concepts  described  in  the  preceding  paragraphs,  the  monitoring  plan  will 
be  implemented,  with  constant  feedback  to  ensure  that  the  model,  the  sampling  plan  and 
the  uncertainty  of  the  monitoring  plan  are  all  considered  when  updating  the  process. 
2.2  Tiered  System  Structure 

Using  the  guidelines  for  establishing  a  monitoring  plan  described  above  leads  to  a 
variety  of  monitoring  plans  based  upon  conditional  aspects  of  the  system  being  studied. 
An  approach  used  by  a  team  of  graduate  students  at  MIT  [9].  provided  a  tiered  system 
approach  to  a  monitoring  plan  that  allows  a  manager  to  best  allocate  the  available 
resources  when  developing  the  monitoring  plan. 


Tier  I 
Indicators 


Tier  II 

Identification  of 
Perturbation 


Tier  III 
Sources  of 
Perturbation 


Water  Column       Living  Resources    Sediment/Benthos 


A/A  1 


DDDDDDD 


Figure  2-  1:  Three  Tiered  Strategy 


The  tiered  structure  is  classified  to  characterize  three  potential  areas  of  impact 
within  a  monitoring  plan  area.  Those  areas  are:  1)  Water  Column  Region,  2)  Living 
Resources  ,  and  3)  Sediment/Benthos.  The  monitoring  plan  addresses  these  three  areas  in 
three  tiered  levels  of  measurement.  The  three  tiered  strategy  is  graphically  represented  in 
figure  4-1.  Each  box  represents  a  variable  that  is  part  of  the  generic  design  study.  The 
number  of  the  boxes  would  change  based  upon  the  specific  application  The  arrows 
represent  that  the  tiers  are  logically  linked  to  each  other.  Note  that  the  logical  links  from 
one  tier  to  the  next  may  indicate  more  than  one  cause  and  effect  relationship  for  a 
particular  variable. 

The  first  tier  in  the  monitoring  structure  is  comprised  of  variables  that  react 
quickly  to  the  changes  in  the  marine  system,  and  provide  advanced  warning  of  potential 
events.  Tier  I  items  include  the  measurement  of  water  quality  variables,  such  as 
temperature,  salinity  and  dissolved  oxygen  in  the  water  column.  The  living  resources 
classification  is  measured  in  fish  population  studies  for  tier  I.  The  benthos  region  tier  I 
measurements  include  bottom  observations  that  may  be  made  with  remote  sensing 
equipment. 

If  the  variables  measured  in  the  first  tier  exceed  a  threshold  value,  then  more 
extensive  analysis  must  be  performed.  This  triggers  the  investigation  of  the  tier  II 
variables  via  a  logical  link  to  the  most  probable  source  of  the  perturbation.  The  tier  II 
variables  require  a  more  complex  analysis  that  is  generally  more  expensive  and  not  be 
performed  in  situ.   An  example  of  a  tier  II  variable  is  the  chemical  analysis  of  the  water 
column  that  requires  a  grab  sample  be  obtained  and  transported  to  a  laboratory.  The 
living  resources  tier  II  variable  would  be  the  analysis  of  animal  or  plant  material  to 
determine  its  make-up  and  compare  it  to  standard  information.   A  tier  II  benthos  region 
variable  would  again  be  a  chemical  analysis  of  the  material  of  the  sediment. 

The  final  portion  of  the  tiered  system  is  the  logical  link  from  the  identification  of 
the  perturbation  to  the  most  probable  source  of  the  change.  This  involves  the  tracing  of 
the  cause  of  the  perturbation  to  its  source.  The  basis  for  this  would  be  set  up  in  the  design 
study,  when  potential  sources  of  perturbation  to  the  marine  system  would  be  identified  as 
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any  facility  that  could  potentially  add  material  and  or  heat  to  the  region  around  the 
boundary  of  the  marine  system  being  monitored  (power  plants,  industrial  plants,  chemical 
plants).  This  detective  effort  is  the  final  resuit  of  the  monitoring  plan  because  it  is  the 
information  that  managers  need  to  form  their  reasoned  opinions  concerning  liability  and 
potential  litigation. 

The  most  significant  advantage  of  the  three  tiered  system  is  its  modularity.  The 
system  can  promptly  be  updated  to  reflect  the  information  obtained  when  the  monitoring 
system  goes  into  effect.  Additionally,  the  three  classifications  of  water  column,  living 
resources  and  benthos/sediment  characterize  the  essential  aspects  of  the  marine 
environment.  Each  tier  can  be  expanded  horizontally  as  information  is  obtained  and  more 
classifications  can  be  added  if  required.  If  additional  links  are  provided  with  the 
improvement  in  knowledge,  the  links  can  also  be  updated. 

The  modularity  of  the  system  allows  the  resource  manager  to  best  use  the  assets 
that  are  available  if  the  entire  plan  can  not  be  fully  implemented.  The  logical  links  that 
connect  the  tiers  describe  flow  paths  that  can  be  chosen  in  either  a  bottom  up  or  top  down 
approach  to  investigate  the  potential  cause  and  effect  relationships. 

2.3  Implementation 

The  use  of  moored  oceanographic  systems  to  measure  components  that  are  part  of 
the  monitoring  plan  will  likely  be  used  because  of  the  cost  effectiveness  of  remote  sensing 
and  the  possibility  of  real  time  data  management  from  a  shore  based  office.  In  order  to 
technically  achieve  the  use  of  moored  oceanographic  systems,  it  is  imperative  that  external 
forcing  functions  that  move  the  system  be  controlled  to  allow  sensitive  equipment  to 
operate  without  damage. 

The  use  of  the  design  process  to  design  a  dynamic  absorber  can  be  used  most 
effectively  in  the  tier  I  analysis  of  the  monitoring  plan.  The  vision  is  that  the  particular 
measurement  desired  to  be  measured  in  the  water  column  in  situ  can  be  designed  to 
operate  installed  on  a  moored  buoy  system.  The  power  source  for  the  system  can  be 
designed  to  last  for  a  designated  period  of  time  that  the  user  would  correspond  to  a 
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maintenance  schedule  for  the  equipment.   Current  renewable  energy  systems  rely  on  solar 
energy  or  other  means  that  would  be  installed  in  the  moored  system  near  the  free  surface 
that  would  suffer  if  excess  drag  or  vibrations  caused  the  buoyant  portion  of  the  moored 
system  to  go  underwater.    Communications  with  the  land  based  researchers  would  also  be 
connected  to  the  installed  devices  via  the  mooring  cable.  This  communication  suite  would 
most  likely  be  installed  in  a  buoyant  body  near  the  free  surface  in  a  vulnerable  position. 

A  specific  example  of  a  typical  oceanographic  system  would  include  acoustic 
transducers  to  monitor  various  components  of  the  marine  environment.  These  transducers 
are  very  sensitive  to  excessive  vibrations,  and  the  design  process  to  remove  as  much 
vibration  energy  from  the  system  as  possible  would  help  to  ensure  that  the  installed 
equipment  could  operate  properly. 

With  the  implementation  of  Autonomous  Underwater  Vehicles  (AUVs)  in  the 
marine  environment,  certain  tier  II  variables  can  also  be  measured.  Ongoing  work  within 
the  Ocean  Engineering  Department  at  MIT  is  investigating  the  implementation  of  new 
technology  that  will  greatly  reduce  the  amount  of  power  required  to  operate  the 
propulsion  systems  of  AUVs.  The  computer  program  required  in  order  to  operate  an 
AUV  in  the  underwater  marine  environment  has  also  been  the  subject  of  recent  studies  at 
MIT.  The  power  and  communication  source  for  the  AUV  would  be  connected  to  another 
moored  cable  system  that  would  have  an  underwater  docking  port.  This  docking  port 
would  provide  the  communications  link  with  the  land  based  researchers  and  replenish  the 
power  used  in  the  tier  II  investigations. 

Of  particular  concern  to  each  of  the  applications  described  above  is  the  principal  of 
flow  induced  vibrations.  Even  with  a  steady  current  moving  past  a  cable  system  or  a  cable 
system  being  towed  by  a  vessel  at  a  constant  speed,  oscillatory  motions  are  induced  when 
vortex  shedding  occurs  on  the  cable.  The  motions  of  the  cable  are  dominated  by  motions 
that  are  transverse  to  the  direction  of  the  flow  and  occur  at  a  frequency  equal  to  the  vortex 
shedding  frequency.  In-line  motions  occur  at  twice  the  vortex  shedding  frequency,  and 
are  not  considered  in  this  analysis  because  of  their  small  contribution  to  the  cable  motion. 
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The  vortex  shedding  excitation  frequency  often  coincides  with  a  resonant  frequency  of  the 
cable  and  results  in  substantial  response. 

The  non-dimensional  parameter  that  characterizes  the  frequency  of  vortex 
shedding  is  the  Strouhal  number,  and  the  typical  value  associated  with  it  is  .21  and  it  is 
constant  over  a  wide  range  of  Reynolds  number.  [4].  The  expression  relating  the  Strouhal 
number  and  the  vortex  shedding  frequency  is: 


St  V 


2n 


d  EQN  2-  1 

Where  co  is  the  angular  frequency,  V  is  the  velocity  of  the  cable  being  towed  or  the  current 
flowing  past  the  cable,  and  d  is  the  diameter  of  the  cable.  In  the  design  process,  a  limit 
needs  to  be  established  for  reasonable  expectations  of  vortex  shedding  frequencies  that  the 
moored  system  may  face.  Assuming  a  limit  of  a  5  knot  current  (2.54  m/sec)  and  a  range 
of  cable  diameters  from  1/4  inch  to  1  inch,  Figure  2-2  represents  the  expected  vortex 
shedding  frequencies. 


0.005 


0.01 


0.015 
Diamtcr  (m) 


0.02 


0025 


Figure  2-  2:   Vortex  Shedding  Excitation  Frequencies  in  a  2.5  m/sec  Current 

In  addition  to  this  vibrational  lock  in  of  the  cable,  an  additional  concern  for  flow 
induced  vibration  is  the  increase  in  the  nondimensional  drag  coefficient,  CD.  This  increase 
in  drag  has  major  implication  to  mooring  design. 


13 


Understanding  the  need  for  the  oceanographic  monitoring  and  the  environment  in 
which  the  components  will  be  operated  leads  directly  to  an  investigation  to  reduce  the 
amount  of  energy  that  the  installed  components  will  see.  This  investigation  to  limit  the 
effect  of  the  ocean  loading,  in  particular  the  vortex-induced  vibrations,  is  pursued  using  an 
in-line  dynamic  absorber. 
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Chapter  3 

Transfer  Matrix  Evaluation 
3.1  Background 

The  transfer  matrix  method  is  a  technique  to  evaluate  a  complicated  system 
through  the  use  of  discrete  transformations  from  one  location  to  another.  For  a  sequential 
system  of  cable,  rigid  bodies,  springs  and  dampers,  the  final  system  transfer  matrix  is  the 
product  of  each  transfer  matrix  from  one  location  to  the  other.  The  benefit  of  this  system 
is  the  ability  to  evaluate  the  response  of  a  complicated  system  through  sequential 
multiplication  of  finite  transfer  matrices  from  point  to  point.  The  analysis  of  this  thesis 
will  evaluate  the  transfer  matrix  for  an  oceanographic  moored  system  with  a  damping 
connector  installed  in-line,  shown  in  Figure  4-1. 

Using  the  transfer  matrix  method  requires  that  a  state  vector  be  defined.  The  state 
vector  is  comprised  of  Fourier  transforms  of  displacements  and  the  corresponding  internal 
forces  [5].  In  the  case  of  a  cable  system  with  installed  masses,  springs  and  dampers,  the 
state  vector  is  described  by: 


Z;    - 


Ty' 


EQN  3-  1 


where  y  is  the  Fourier  transform  of  the  displacement  of  the  system,  y'  is  the  Fourier 
transform  of  the  rotation  of  the  system,  and  T  is  the  tension  of  the  system.   All  three  items 
are  evaluated  at  the  particular  point  (annotated  by  the  subscript  /),  and  the  transformation 
from  one  location  to  the  next  is  defined  by: 

z,  =  4.2,_,  EQN  3-  2 

where  the  transfer  matrix  A.  links  the  two  state  vectors. 
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3.2  Derivation  of  a  Transfer  Matrix 

The  determination  of  the  transfer  matrix  is  based  upon  an  analytical  description  of 
the  cable  system.  This  thesis  uses  the  transfer  matrix  library  determined  by  Li.[l].  In 
order  to  describe  how  the  transfer  matrices  are  established,  an  alternate  description  of  an 
equivalent  system  is  provided.  Knowing  the  geometric  relationships  and  determining  the 
equations  of  motion  are  the  keys  to  deriving  the  transfer  matrix  for  a  system. 


Figure  3-1:    Schematic  of  System 

Using  figure  3-1  and  setting  Mi=M?-M,  the  following  equations  of  motion  can  be 
described  in  the  y  direction: 


M(yc2+ycl)=T1/-Tl_ly;_ 


EQN  3-  3 


Where  T    is  the  internal  force  at  the  joint  located  at  y0  acting  in  the  yo  direction 

The  equations  for  summation  of  moments  about  the  center  of  gravity  of  the  two 
links  are  described  by  the  following: 
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Ij2  =  A'„(0,  -  02)  +  KiO,  -$2)+  T,(y;  -  02)f  EQN  3-  4 

I  J,  =  Ke{62  -$,)+  R802  -6,)  +  77_10;_]  -  9x)d  EQN  3-  5 

Where  K$  is  the  value  of  the  torsional  spring  constant  (Nm/rad)  and  Rq  is  the  value  of  the 
torsional  damping  constant  (Nm/(raaVsec)). 
If  the  following  geometric  relations: 

yc>=yy,->+yyo  eqn3-6 

Ly  Ly 

V«  =  J-yQ  +-J-*  EQN  3- 7 

are  inserted  in  Equation  3-3,  we  obtain: 


?*  =  75  ^a/^-^'^ZTT;^-'  +T^  eqn 3- 8 

(d  +  f)M  (d  +  f)   U  Li 


And  if  we  substitute  the  following  angular  relationships: 

Q  =  hZllzL    and       ^    =  ZlZA  eqN  3.  9 


into  equations  3-4  and  3-5,  we  obtain: 

EQN  3-  10 

and 

EQN  3- 11 
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Performing  the  Fourier  transform  of  Equation  3-8  and  setting  Ti=Tj.i=T,  d=f=l2 ,  c=e=li, 
and  Li=L2=L,  we  obtain: 

7  =  — tttP-  -7,-,]-^,-,  +7)  EQN  3-  12 

a)    2/:  l  J     2/2 

Where  y0,  yt_x ,  yi ,  >'',_, ,  and  >>',    are  the  Fourier  transforms  of  the  displacements  and 
the  slopes  at  the  particular  locations. 

Performing  the  Fourier  transform  of  Equation  3-10  and  Ic2  =  Mk2  (k  is  the  radius 
of  gyration),  we  obtain:    • 


J/rW)-^  -y.)  =  A-9^o-(.v„,+y,)  +  ^M2y0-Cv<„+y,)1  +  r->  _y±-y±  ^       EQN  3.  13 

Performing  the  Fourier  transform  of  Equation  3-11,  and  setting  Ici  -  Mk2  (k  is  the  radius 
of  gyration),  we  obtain: 


14 

If  we  substitute  Equation  3-12  into  Equation  3-13  and  rewrite,  we  obtain: 

Adi  +5,7.  =Ol.,+Dj,_1  EQN  3- 15 

Where: 

4  = (1  +  — )  +  — -  +  — -  +  - — -  +  '   g       +  -=—  +  J— 

'  L  2l2        L       LU         L  Ll2  L      1L 


2         1..2 


2/2         /,&r     yo>/2  2<y 

MkW_       KL_l]K0     jcoR9     l^Rgja) 


2Z/2  L       Ll2         L  Ll2 
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a  = 


Mk-T     TKg      R^T      77,        T2 


2/,        /..a;2      7<y/2      2L     2co2 


And  combining  Equations  3-13  and  3-14  and  rewriting  to  obtain: 


A2y,  +Bzy,  =Ciy,~\  +  Ajv 


EQN  3-  16 


Where 

-Mk2co2      Tl2  -Mk2co2      Tl2 

A,  = +  — -    and    B-,  =  -  TL    and    C,  = and  Z),  =  77, 

L  L  2  2  L  L 

The  final  equations  for  the  transfer  function  can  be  formed  using  the  following 
relationship: 


y± 


Ml        M2 


yj-i 


Thus: 


'n  = 


2?2C,  —  5,C; 
^  -A2B, 


'»  = 


1  B2D,-BXD2 

T  A,B.  -  A,B, 


f  ™]^2  2^1 


/4,52  -/i:5, 


22    - 


/J,Z)2  -  /J:£>, 
A,B2-  A2B} 


EQN  3-  17 


Using  the  methods  of  transfer  matrix  operations  allow  an  efficient  calculation  of 
important  parameters  such  as  the  amount  of  energy  that  is  reflected,  absorbed  or 
transmitted  through  a  system.  The  analysis  in  the  next  chapter  will  provide  that 
information  and  the  techniques  needed  in  order  to  reduce  the  amount  of  flow-induced 
vibrations  that  are  seen  in  the  system 
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Chapter  4 


In-line  Dynamic  Absorber 

4.1  System  Description 

Using  the  principles  described  in  chapter  3,  the  transfer  matrix  formulation  for  a 
specific  system  can  be  described  and  studied.  The  system  studied  is  a  hinged  connection 
of  two  rigid  lump  masses  connected  with  a  linear  translational  spring  and  damper 
connected  in  an  infinitely  long  cable,  shown  in  figure  4-1.  The  cable  is  aligned 
horizontally  and  the  tensions  are  assumed  to  be  constant  along  the  cable. 


M2,I 


2,-12 


h=c+d 
l=e+f 


c     d       , 

c           c            K          e            .  .       f       . 

1 

Figure  4-  1:  Diagram  of  System 

Where  the  relationship  between  the  torsional  components  defined  in  chapter  3  and  the 
linear  translational  components  described  here  are  given  by  the  following: 


Kg  =  Ka1        and        R^  =  Rb: 


EQN  4-  1 


An  incident  wave  traveling  along  the  cable  can  be  described  by: 


-  jkx     jqM 


Ae-Jlae 


EQN  4-  2 


where  A  is  the  amplitude,  k  the  wave  number  and  co   the  radial  frequency  of  the  wave. 
The  subscript  /  corresponds  to  the  incident  wave. 

When  the  incident  wave  encounters  a  discontinuity  along  the  cable,  the  wave  can 
be  reflected  or  transmitted  through  the  region.   A  comparison  of  the  wave  amplitudes 
provides  the  following: 
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71  =  -*-    and    R=^-  EQN  4-  3 

4  A 


where  7^  is  the  transmission  coefficient  and  Rc  is  the  reflection  coefficient. 

Knowing  the  rate  at  which  energy  is  carried  by  an  incident  wave  along  the  cable  is 
defined  by 

E   =-pcCoj2A  2  EQN  4- 4 

mm 

Where  pc  is  the  linear  density  of  the  cable,  C  is  the  wave  speed,  co  is  the  frequency  and  A 
is  the  amplitude  of  the  wave. 

When  the  wave  encounters  the  system  shown  in  figure  4-1,  the  energy  of  the  wave 
can  be  reflected,  transmitted,  or  absorbed  into  the  damper  between  the  two  masses.  The 
energy  of  the  waves  that  reflect  back  from  the  system  can  be  determined  as  well  as  the 
energy  of  the  waves  that  transmit  through  the  system.  The  remaining  energy,  if  any, 
represents  the  energy  absorbed  by  the  system,  considering  that  energy  must  be  conserved. 
The  ratio  of  these  energies  are  described  by 


E. 


E.. 


A 


P=—   alternatively  p  =    — H  EQN  4- 5 


: 


KAJ 


and 


E,     ..         .    ,  Z, 


f 


E,  Z, 


A, 


y  =  — '-  alternatively  y  =  —  —  EQN  4-  6 


A,  J 


therefore 


£ 

a  =  —  alternatively  a  =  \- p-y  EQN  4-  7 

E, 


where  the  subscript  /  stands  for  incident,  r  stands  for  reflected,  /  stands  for  transmitted, 
and  a  stands  for  absorbed. 

But  the  reflected  and  transmitted  power  ratios  can  also  be  described  in  terms  of  the 
reflection  and  transmission  coefficients  (equation  4-3). 
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p  =  \Rc\-  and  r  =  ^-|re|2  EQN4-8 


where  a  is  still  described  as  in  equation  4-7 


4.2  Resonance  Frequencies 

The  first  step  of  the  analysis  is  to  determine  the  frequency  at  which  the  system 
would  have  the  maximum  amount  of  energy  absorption.  This  occurs  at  a  frequency  when 
the  reflection  coefficient  is  zero,  thus  all  of  the  energy  is  either  transmitted  or  absorbed. 
Using  the  elements  of  the  transfer  matrix,  the  expression  for  the  reflection  coefficient  is 
provided  by: 

-    (t2i(<»  ))+(<■>)    Z  ,  Z2t  ]2(w)  h-j  -w(z  j  t  22(cd  )-  Z  -,  t  n((J)< 
Rc(w)   =— 


t21(0))-(u>)    Z  !  Z  2  t  12(o>)+j    <a-(z  ^t^uJ  +  Z  2t  u(»))  EON  4-  9 

Where  the  elements  of  the  square  transfer  matrix  are  tij(co),  for  i  and  j  equal  to  1  to  2  and 
the  impedance  of  the  input  and  output  cable  are  Z\  and  Z2,  respectively. 

In  addition  to  the  reflection  coefficient,  the  transmission  coefficient  can  also  be 
expressed  in  terms  of  the  transfer  matrix.  The  expression  for  the  transmission  coefficient  is 
provided  by: 

j  -2-ta-Z  , 
Tc(o>)   = 


2 


Rc(u)     = 


t  2i(ft))-((o)    Z  t  Z  2  t  12(  (■))•+- j    <d(Z  ,  t22((D)-t-Z2  t  u(0)))  EON  4-10 

Reducing  the  expression  for  the  reflection  coefficient  aids  in  the  location  its  zeros: 

f  2422  222  222  2  32  ^ 

(-t  2](w)    +0)    Z  ,    Z  j    t  j2(0))   •+-(■)    Z  j    t  ^-Ju>)    -  0)    Z  2    t  |  |((i>)    +  2i  <o  Z  j  t  ,^(  u  )  t  ->](<i> >  +  2i  u    Z  ,  Z  2    t  |  |(o>)  t  i2(<*)/ 

/  2  2  42^  "'22  2  ""  2       2  -\ 

it  2j((D)   -  2  t  2](u)  0)    Z  j  Z2  t  |2Cfl)-f-w    Z  |    Z  2~  t  12((i>)   +  o>    Z  ]    t  2t(w)    +~  -  u  'z  |  '  22^  "  '  Z  2  '  1 1(  w  '  +  u    z  2    l  ll'0"'  / 

EQN4-11 
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The  numerator  of  this  expression  can  be  solved  for  values  of  co  such  that  it  equals  zero. 
These  values  of  co  are  the  frequencies  at  which  there  is  no  reflection  and  the  system  has  the 
potential  for  maximum  energy  absorption. 

A  second  method  of  calculating  the  frequencies  of  the  zeros  of  the  reflection 
involves  the  use  of  the  input  impedance,  defined  by: 


_        force 


velocity 


and      z_(a) 


t21(o))-hj    «aZ2t  n(w) 

2„ 
(i)  -Z  2t  12(w  )H-j  -<o-t  22^ u  > 


EQN  4-  12 


By  using  the  principal  of  impedance  matching,  the  values  of  co  such  that  the 
imaginary  part  of  this  impedance  is  zero  will  provide  the  values  of  the  frequencies  where 
the  potential  for  maximum  energy  absorption  will  occur.  This  leads  to: 


Imaginary  part  of    z  m(<»)  = 


-a   -Z2    t  ij(u)-t  12((u)-t  21((D)-wt22(u) 

4,2  2        2  2 

a>    Z2    t  12((i))   -(-a)    t  22((i) ) 


EQN  4-  13 

The  numerator  of  this  expression  can  also  be  solved  in  order  to  find  the  values  of  co  such 
that  the  expression  equals  zero.  These  values  are  the  resonance  frequencies  of  the  system. 


Table  4-  1:  Input  Values  for  the  Determination  of  the  Reflection  Coefficient  Zeros 


Characteristi 
c 


M, 


I. 


2l 


TV, 


K 


H 


Value 


.0254  m 


1524  m 


.1778  m 


4.52  kg 


.035 


1.144  kg/m 


1000  newton 


10  newton/m 


1000  newton 


— ■ 


Characteristi 
c 


M, 


£i 


T, 


R 


Value 


.0254  m 


1524  m 


.1778  m 


4.52  kg 


.035  kgm* 


1.144  kg/m 


1000  newton 


0  newton/(m/sec) 


An  analytic  approach  and  a  corresponding  numerical  solution  to  this  problem  is 
outlined  in  Appendix  A.  With  the  value  of  the  damper  set  to  zero,  a  numerical  solution 
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can  be  obtained  for  both  the  zeros  of  the  input  impedance  expression  and  the  reflection 
coefficient  expression.  A  numerical  solution  to  this  problem  requires  the  input  of 
characteristic  values  of  the  physical  parameters  of  the  system.  In  order  to  compare  the 
simulation  with  a  real  world  model  of  a  cable  system  modeled  by  Vandiver  [6],  the  values 
of  the  system  shown  in  Figure  4-1  were  set  in  accordance  with  Table  4-1 . 

The  specific  values  shown  in  Table  4-1  represent  a  real  world  example  of  some 
oceanographic  mooring  components.  Using  these  characteristic  properties,  a  computer 
simulation  of  the  system  was  conducted  and  the  value  of  the  frequencies  of  the  zeros  of 
the  reflection  coefficient  were  obtained  using  the  secant  method.  The  resulting  values 
were  Q.\  =  42.6  rad/sec  and  Q2 =  94. 1  rad/sec.  Using  the  impedance  method,  the  same 
values  were  obtained.  The  values  of  the  reflection  coefficient  and  the  transmission 
coefficient  are  plotted  in  Figure  4-2. 


250 


(radians  sec) 


Figure  4-  2:  Graph  of  Reflection  and  Transmission  Coefficients  Using  Input 

Parameters  from  Table  4-1  (R=0) 

For  the  curve  shown  in  figure  4-2,  the  value  of  the  damper  was  set  at  zero  to 
observe  the  location  of  the  minimum  values  of  the  reflection  coefficient  with  no  energy 
being  absorbed  in  the  system  to  aid  in  determining  a  starting  value  for  the  secant  method 
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analysis.  The  location  of  these  minimum  values  is  not  affected  by  a  change  in  the  damper 
value,  R,  although  the  values  of  the  reflection  and  transmission  coefficients  are.  This  is 
shown  in  Figure  4-3,  by  increasing  the  damper  to  500  N/(m/sec).   The  difference 
represents  the  amount  of  energy  that  is  absorbed  in  the  system 
1 


(radians*  sec) 


Figure  4-  3:  Graph  of  Reflection  and  Transmission  Coefficients  Using  Input 
Parameters  from  Table  4-1  (R=500  N/(m/sec)) 


4.3       Optimum  Value  for  the  Damper 

The  determination  of  the  optimum  value  of  a  damper  to  absorb  the  maximum  value 
of  energy  is  frequency  dependent.  The  maximum  amount  of  energy  will  be  absorbed  at 
the  resonance  frequencies  determined  in  the  previous  section.  Using  the  two  resonance 
frequencies  as  a  basis,  a  search  was  conducted  to  determine  this  optimum  damper  value 
using  the  same  expressions  outlined  in  Appendix  B.  This  time,  however,  all  of  the 
equations  are  dependent  upon  the  value  of  a  particular  damper  with  the  frequency  held 
constant  at  either  Qi  or  Q2    This  investigation  emphasizes  the  higher  resonance  frequency 
because  there  is  more  relative  motion  between  the  two  masses  that  leads  to  more  energy 
absorption  in  the  damper  as  seen  in  the  modal  sketch  in  Appendix  A. 


25 


In  order  to  determine  the  value  of  the  damper  that  will  lead  to  the  maximum 
amount  of  energy  absorbed  for  a  particular  frequency,  the  absorption  power  ratio  must  be 
analyzed.  Using  equation  4-6,  the  derivative  of  the  absorption  power  coefficient  with 
respect  to  R  is  described  by: 

5R  5R     5R  EQN  4-  14 

Substituting  equation  4-7  into  equation  4-13,  the  following  is  obtained: 

1_ 
8R  5R  6R  EQN  4-  15 

Which  can  be  rewritten  as: 


5a     _     -HRcRCj        \Z1      C     7 


— \      8 T    -T 


5a  5Rc—  5Rc       z2  f6Tc-  5T 


5r  Wr     c      c  5r  /    zi  Wr     c      c  5r  /  EQN  4- 16 


Because  the  order  of  the  partial  derivative  operation  and  the  complex  conjugate 
calculation  are  interchangeable,  the  final  expression  can  be  rewritten  as: 


fa                8R  =  ^     o      M        Z2    5Tc—     T      /5TcVl 
—     - Rc+Rc Tc  +  Tcl 

Sr  \5r      c        \5r//    zi  \5r     c         \5r//  EQN  4- 17 

The  partial  derivative  of  the  transmission  coefficient  and  the  reflection  coefficient 
require  that  the  partial  derivative  of  the  transfer  matrix  with  respect  to  R  be  obtained. 
This  analysis  is  described  in  Appendix  B  and  is  based  upon  the  work  of  Li  [  1  ].  The 
analysis  provides  expressions  for  the  two  partial  derivatives: 

(02Z:(ll:(R)t:lj^e(R)-t|2^rme(R)t:i(R)  +  tn_pnme(R)t::(R)-l;lkR)t;:j)nme(R)) 

+j  {<o32::(t1,_pnme(R)t,2(R)-tn(R)t12j)rime(R))-a)  (t  21jillll{R)-i  2:(R)-  t  :](R)t  22_prime(R)) 


5RC    _ 


5R  r/  2  \  ]2 

^t2,(R)-u*Z  ,  Z2i  ,2(R)j+j    w  (t22(R)Z  ,  ri  ,^R)  Z  ,)j 


EQN  4-  18 
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and: 


5R 


2Z, 


tt)2-(t22j^mctR)Z1^tnj?nn,e(R)Z2)+j  •  (M3-Z  ,Z  2  t  l2jtimt[  R)  -  »  t  21j>finte(R)) 
t21(R)-a)    Z  r-Z  2t  i2<R)j-t-J    0)^t22(R)Z  ,-H  u(R)Z2)j 

EQN  4-  19 


The  determination  of  the  optimum  value  of  the  damper  to  absorb  the  maximum 
energy  becomes  a  numerical  solution,  given  the  physical  characteristics  of  the  mooring 
system.  The  extreme  case  is  the  replacement  of  the  second  mass  to  encounter  the 
incoming  wave  in  Figure  4-l(M2)  with  a  hinge  termination.  The  optimum  value  of  the 
damper  would  be  when  the  damper  impedance  was  equal  to  the  impedance  of  the  input 
cable: 

11  R  Rb2 

or  Damper  Ratio  = 

^■optimum  P^ 


Ropamim  =  pc—  therefore:    Damper  Ratio  = or  Damper  Ratio  = 


EQN  4-  20 

Where  R  is  the  value  of  the  damper,  b  is  the  transverse  distance  from  the  centerline  of  the 
system  to  the  location  of  the  damper,  L  is  the  length  of  the  first  mass  (  h  in  Figure  4-1  ),  p 
is  the  linear  density  of  the  cable,  and  c  is  the  speed  of  the  wave  traveling  along  the  cable. 
A  comparison  with  the  value  of  this  optimum  damper  will  provide  a  great  deal  of  insight 
on  the  performance  of  the  evaluation  and  design  process. 

The  design  process  described  in  Appendix  B  was  based  upon  the  fact  that  the  user 
will  know  the  value  of  the  linear  density  of  the  cable  and  the  range  of  tension  expected  in 
the  system.  From  these  two  parameters,  the  impedance  of  the  cables  (Zi  and  Z2  in  Figure 
4-1)  can  be  easily  determined.  Using  those  impedances  as  the  entry  point  in  the  design 
process,  a  new  set  of  resonance  frequencies  must  be  obtained  using  the  process  described 
in  Appendix  A.  Based  upon  the  highest  of  the  two  frequencies  calculated,  the  optimum 
value  of  the  damper  is  numerically  determined  using  the  secant  method  for  the  derivative 
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by  substituting  equations  4-19  and  4-18  into  equation  4-17.  This  process  is  done  over  a 
range  of  values  for  the  cable  impedance,  and  the  results  are 
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Figure  4-  4:  Plot  of  Optimum  Damper  Value  versus  Cable  Impedance  Using 
Input  Parameters  from  Table  4-1  and  Evaluated  at  the  Corresponding  Cl2 

plotted  in  Figure  4-4.  The  range  of  the  cable  impedance  varied  from  10  to  200  kg/sec, 
maintaining  the  value  of  the  linear  density  of  the  cable  constant.  The  linear  density  as  well 
as  the  rest  of  the  physical  characteristics  are  the  same  as  the  values  as  outlined  in  Table  4- 
1,  with  the  exception  of  the  tension  at  T,  and  T;.i  which  must  be  altered  in  order  to 
maintain  the  desired  value  for  the  cable  impedance  maintaining  the  linear  density  of  the 
cable  constant. 

A  plot  of  the  values  of  the  absorption  power  ratio  versus  the  damper  values  for 
varying  increments  of  the  cable  impedance  is  provided  in  Figure  4-5.  When  the  damper 
values  are  divided  by  the  cable  impedance  for  the  particular  variant  (Zi=Z2=Zc),  the 
curves  closely  follow  each  other,  with  a  maximum  absorption  power  ratio  of 
approximately  .5  at  a  damper  ratio  of  approximately  .11.  This  is  shown  in  Figure  4-6. 
The  variance  of  the  cable  impedance  does  not  have  a  drastic  affect  on  the  maximum 
absorption  power  ratio. 

The  variation  of  the  optimum  damper  is  also  affected  by  the  mass  of  the 
component  masses  Mi  and  M2.  With  the  other  physical  characteristics  held  constant,  the 
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Figure  4-  5:   Plot  of  Absorption  Power  Ratio  for  Varying  Input  Impedances  Using 
Input  Parameters  from  Table  4-1  and  Evaluated  at  the  Corresponding  Q2 
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Figure  4-  6:  Plot  of  Absorption  Power  Ratio  Using  Input  Parameters  from  Table  4-1 
and  Evaluated  at  the  Corresponding  il2  for  Varying  Input  Impedances 
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Figure  4-  7:  Plot  of  Absorption  Power  Ratio  for  Varying  Mass  Ratios  (Mi=M2) 
Using  Input  Parameters  from  Table  4-1  and  Evaluated  at  the  Corresponding  Q.2 


change  in  the  optimum  value  of  the  damper  with  the  varying  values  of  the  masses(Mi=M2) 
is  shown  in  Figure  4-7.  Here  the  mass  of  the  two  components  is  compared  to  the  mass  of 
an  equivalent  length  of  cable.  The  mass  ratio  is  defined  by: 


MR  = 


A/,  +M2 
2P]d  +  2p2f 


EQN  4-  21 


Note  that  the  increase  in  mass  tends  to  increase  the  absorption  power  ratio,  and  it 
asymptotically  approaches  a  maximum  value  of  .5  when  both  masses  are  increased  the 
same  amount.  If  the  two  masses  are  not  altered  equally,  varying  the  mass  ratio  will  alter 
the  absorption  power  ratio  differently.  Figure  4-8  shows  that  an  increase  in  the  mass  of 
Mi  while  maintaining  M2  constant  tends  to  decrease  the  value  of  the  absorption  power 
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Figure  4-  8:   Plot  of  Absorption  Power  Ratio  Varying  Mass  Ratios  (Alter  Mi  only) 
Using  Input  Parameters  from  Table  4-1  and  Evaluated  at  the  Corresponding  Cl2 
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Figure  4-  9:  Plot  of  Absorption  Power  Ratio  Varying  Mass  Ratios  (Alter  M2  only) 
Using  Input  Parameters  from  Table  4-1  and  Evaluated  at  the  Corresponding  Q2 
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Figure  4-  10:   Plot  of  Absorption  Power  Ratio  varying  Mass  Ratio  (Alter  Mi  or  M2onIy) 
Using  Input  Parameters  from  Table  4-1  and  Evaluated  at  the  Corresponding  Q2 

ratio,  while  Figure  4-9  shows  that  an  increase  in  the  mass  of  M2  while  maintaining  Mi 
constant  tends  to  increase  the  value  of  the  absorption  power  ratio.  Figure  4-10  shows  the 
trend  towards  reaching  the  ultimate  value  of  total  absorption  at  the  second  resonance 
frequency  as  the  mass  ratio  of  the  system  is  increased  by  altering  the  mass  ratio  with 
increasing  Mi  or  M2  only. 

Another  physical  characteristic  of  the  system  involves  the  lengths  of  the 
components.  If  the  lengths  of  the  components  are  changed,  the  natural  frequencies  and 
the  corresponding  values  of  the  optimum  damper  will  also  change.  Figure  4-1 1  represents 
the  effect  that  changing  the  length  of  Mi  and  M2  has  on  the  performance  of  the  system, 
using  the  following  relationship  for  the  length  ratio: 


LR  = 


f 


EQN  4-  22 


Where  d  and  fare  the  lengths  described  in  Figure  4-1 .  Note  that  the  location  of  the  center 
peak  for  the  value  of  the  damper  in  Figure  4-11  is  not  constant  for  the  five  different  length 
variations  as  it  was  for  the  mass  and  the  impedance  variants.  An  increase  in  the  length  of 
Mi  while  keeping  the  length  of  M2  constant  causes  a  decrease  in  the  absorption  power 
ratio,  whereas  the  increase  in  length  of  M2  while  keeping  Mi  constant  increases  the 
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Figure  4-11:   Plot  of  Absorption  Power  Ratio  Varying  Length  Ratio  Using  Input 
Parameters  from  Table  4-1  and  Evaluated  at  the  Corresponding  H2 

absorption  power  ratio.  The  geometry  of  the  system  is  the  dominant  term  in  determining 
the  location  of  this  peak  value.  While  altering  the  length,  the  other  constants  outlined  in 
Table  4-1  remained  constant  with  the  exception  of  the  moment  of  inertia  which  is  affected 
by  the  geometrical  changes. 

The  common  factor  in  both  the  mass  and  the  length  analysis  is  that  an  increase  in 
the  moment  of  inertia  in  the  first  mass  to  see  the  incoming  wave  as  defined  by  the  analysis 
(Mi),  whether  caused  by  an  increase  of  mass  or  a  geometrical  change,  decreases  the 
absorption  power  ratio,  and  less  energy  is  absorbed  in  the  system.  An  increase  in  the 
moment  of  inertia  of  the  second  mass  in  the  system  to  see  the  wave  (M2),  will 
correspondingly  increase  the  absorption  power  ratio  and  improve  the  system  performance 
if  the  location  of  the  external  force  can  be  identified. 

Even  though  the  lower  resonance  frequency  provides  a  location  when  the 
reflection  coefficient  is  zero,  the  system  does  not  have  enough  relative  motion  between  the 
two  masses  to  absorb  a  great  deal  of  energy.  This  is  shown  in  Figure  4-12. 
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Figure  4-  12:  Comparison  of  Absorption  Power  Ratios  for  the  Two  Resonance 
Frequencies  Using  Input  Parameters  from  Table  4-1 

The  complexities  of  the  algebra  involved  with  the  formulation  of  the  design  tools 
could  lead  to  an  error  that  is  not  readily  noticeable.  An  attempt  was  made  to  validate  the 
results  of  this  analysis  with  an  independent  source.  Hayner  [7]  used  a  spectral  finite 
element  method  to  calculate  transmission  power,  absorption,  and  reflection  in  cable 
systems.  Similar  to  conventional  finite  element  methods,  a  stiffness  matrix  is  formulated, 
which  relates  the  element  displacements  to  the  applied  forces  at  nodes  due  to  either  other 
elements  or  external  forces.  The  stiffness  matrix  includes  mass  and  dissipative  effects,  as 
well,  by  letting  the  matrix  elements  be  complex.  Formulating  the  cable  dynamics  in  terms 
of  a  stiffness  matrix  brings  to  bear  well  known  finite  element  solution  procedures 
discussed  by  Bathe  [8]  for  solving  the  response  of  interconnected  elements.  The  spectral 
finite  element  differs  from  a  conventional  finite  element  in  that  the  formulation  is  based  on 
a  wave  solution,  which  is  exact  to  within  the  assumptions  inherent  in  the  governing 
differential  equation.  Two  additional  matrices  relating  element  displacements  to  wave 
amplitudes  and  forces  to  wave  amplitudes  are  used  to  calculate  reflection,  transmission, 
and  absorption  coefficients.  The  spectral  finite  element  approach  may  also  be  considered  a 
generalization  of  the  transfer  matrix  approach  discussed  by  Fertis  [9]. 
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Figure  4-  13:  Error  Comparison  with  Hayner  Method  of  Transmission  Power  Ratio 
Using  Input  Parameters  from  Table  4-1  and  R=  20  kg/sec  and  Ii=I2=.005  kgm2 

A  comparison  of  the  two  methods  is  provided  in  Figure  4-13  and  4-14.  The  input 
data  is  provided  in  Table  4-1  with  the  exception  that  the  value  for  the  damper  is  set  at  20 
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Figure  4-  14:  Error  Comparison  with  Hayner  Method  of  Absorption  Power  Ratio 
Using  Input  Parameters  from  Table  4-1  and  R=  20  kg/sec  and  Ii=I2=.005  kgm2 
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kg/sec  and  the  moment  of  inertia  of  the  two  masses  ( I  1=  I2)  is  .005  kg  m2.  Figure  4-13  is 
the  absolute  value  of  the  difference  in  the  transmission  power  ratios  for  the  two  methods. 
Figure  4-14  is  the  absolute  value  of  the  difference  in  the  absorption  power  ratio  between 
the  two  methods.  Both  comparisons  reveal  that  the  two  methods  provide  nearly  identical 
results  for  both  the  absorption  power  ratio  and  the  transmission  power  ratio 
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Chapter  5 

Final  Results/  Recommendations  for  Further  Research 

The  analysis  of  the  moored  cable  system  with  the  in-line  absorber  was  very 
effective.  For  varying  values  of  cable  impedance,  an  optimum  value  of  the  damper  was 
determined.  Throughout  the  entire  range  of  the  possible  forced  vibrations  considered 
from  vortex  induced  forces,  there  was  a  corresponding  value  of  a  damper  that  could 
mitigate  the  effect  of  the  motions.  This  will  allow  the  tuning  of  the  system  to  help  prevent 
the  lock-in  of  the  response  of  the  cable  with  a  forcing  function  near  its  resonant  frequency. 

Throughout  the  entire  effort  to  determine  the  value  of  the  damper  that  would 
provide  the  maximum  amount  of  power  absorption,  the  magnitude  of  the  absorption 
power  ratio  was  limited  to  .5  when  alterations  were  made  symmetrically  with  respect  to 
the  characteristics  of  the  two  masses.  This  matches  well  with  previous  attempts  to 
maximize  the  energy  absorbed  in  a  cable  system  by  Vandiver  and  Li  [10],  although  that 
analysis  involved  the  use  of  a  different  type  of  absorber  that  had  three  masses  instead  of 
two  and  two  springs  and  two  dampers. 

When  an  investigation  was  performed  on  the  behavior  of  the  system  with  increased 
cable  impedance  while  maintaining  the  linear  density  constant  the  results  indicated  a  near 
linear  relationship  between  the  value  of  the  damper  required  in  order  to  absorb  the 
maximum  amount  of  energy  and  the  tension  seen  on  the  cable.   As  the  tension  increased, 
the  value  of  the  damper  required  in  order  to  absorb  the  maximum  amount  of  power  also 
increased.  A  second  independent  researcher,  Hayner,  was  able  to  duplicate  the  results 
obtained  in  chapter  four  with  a  basis  of  only  the  input  parameters,  an  indication  that  the 
process  was  correct. 

The  alterations  of  the  physical  characteristics  of  the  masses  had  a  great  affect  on 
the  behavior  of  the  system  as  a  whole.  An  increase  in  mass  of  the  system,  if  divided 
equally  between  the  two  component  masses,  did  not  provide  as  much  absorption  as  a 
concentration  of  the  increase  in  mass  in  M2  alone,  where  an  absorption  power  ratio  of 
nearly  1.0  was  obtained  when  the  mass  of  M2  was  16  times  Mi.  The  limit  to  that  approach 
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is  that  the  differences  between  the  two  masses  would  have  other  second  order  affects  on 
the  system  that  are  not  modeled  in  this  analysis.  An  increase  in  the  mass  of  Mi  alone 
decreased  the  value  of  the  absorption  power  ratio,  decreasing  optimum  power  absorption. 

The  alteration  of  the  component  mass  lengths  also  had  a  significant  effect  on  the 
maximum  power  absorbed.  An  increase  in  the  length  of  Mi  decreased  the  absorption 
power  ratio,  and  an  increase  in  the  length  of  M2  increased  the  absorption  power  ratio. 
These  trends  follow  the  observations  of  the  mass  changes  described  above.  In  both 
instances,  the  changes  reflected  changes  in  the  moments  of  inertia  of  the  component 
masses.  When  the  moment  of  inertia  of  M2  was  increased  such  that  it  resembled  a  fixed 
pin,  the  power  absorption  was  increased. 

The  first  area  that  requires  further  research  is  the  application  of  this  process  to  the 
finite  case,  where  the  existence  of  the  boundary  conditions  at  the  termination  of  the  system 
would  alter  the  results  of  this  infinite  cable  application. 

The  dynamic  response  of  the  system  to  varying  currents  or  changes  in  speed 
should  be  included  in  a  further  study  of  the  system. 

The  inclusion  of  the  bending  stiffness  of  a  cable  would  provide  a  more  robust 
model  with  a  more  detailed  wave  propagation. 

The  construction  and  testing  of  a  real  world  absorber  would  apply  the  principles  of 
this  thesis  and  validate  the  power  absorption  methodology. 

Additionally,  an  investigation  should  be  made  into  the  effects  that  multiple  in-line 
damper  systems  will  have  on  the  behavior  of  the  system  as  a  whole. 
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Appendix  A 
Determination  of  the  Zeros  of  the  Reflection  Coefficient 


R 


z,         Tm 

i 

M,  I, 

B 

If 

Ti 

*-'l 

n —  "  * 

i  ivi2,i2 

i 
J 

h=c+d 

<     d 

>  -           c 

_K_ 

i     ^           f 

l=e+f 

> 

■v 

1 

Figure  A-l:  Schematic  of  the  System 

Based  upon  the  formulations  of  the  transfer  matrix  equations  from  Li,  the 

following  analysis  describes  the  approach  to  determining  the  location  of  the  frequencies 

where  the  power  absorption  will  be  a  maximum.  A  computer  program,  Mathcad  5.0®, 

was  used  in  the  calculation  process. 

The  first  expressions  to  determine  the  transfer  matrix  are  expressed  in  terms  of  the 
values  of  the  particular  dimensions: 

a  j(0))  :=T  j  minut  d+Hc+Ka  -(o>)    I  j  -t-j  •ffl-R-b  a2(w)  :=T  -f+H-e+K-a2- (u )  - 1  -^  -t- j  -es-R-b2 

2  2 

a  3(a>)   -Ki  -)-j    (fl-Rb 

The  Intermediate  equations  are  expressed: 
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The  intermediate  expressions  continue: 
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(a  2(u  )h  +  a  3(01  )-lj 


M  j-d-1 


1  +  Myf2-: 


Now  icnowing  these,  the  transfer  matrix  can  be  determined: 


T  4(«>)-l-T  !(o)-Y  6(®)  Y  5(<")+7  2(»)7  6( «  ) 

l  ll(»)  — T  i(a>)+7  3(0) ; 7- 7- t  12C*)    -7  2(w)-t-T  3(w)- 


1-7  3(w)7  6("0 


1-7  3(0)7  6(») 


■2l(a) 


7  4(<»)+7  !(w)-7  e^) 
1-7  3(0  )7  6(w) 


22(o>) 


7  5(»)+7  2(<l))'7  e^) 
1-7  3(0  )7  6(«) 


Now  the  reflection  and  transmission  coefficients  are  determined 


R  .(«•). =- 


(t2i(  »))+(<»)    Z  i'Z  2*  12<a))j+J  'ffl(Z  it22(w)-Z2t  ,,((.))) 


t21(0))-(a))    Z  ,  Z2t  !2(a))+j    w(Z  ,  t  22(&))  +  Z2t  u(a))j 


Tc(»)  := 


j  -2(oZ  j 

2 

t21(o))-((o)    Z  ,  Z2t  12((i))4-j    <o(z  ,  t22(a))-f-Z2t  n(w)) 


The  Power  Reflection,  Transmission  and  Absorption  Ratios  are: 


P(»):=R   (u)-R   («o)  y(<o)  := — Tc((i))Tc(u)  a(w)  .  =  1-P(m)-y(4)) 


1 


The  input  impedance  is: 


„      -    force  .  ,  t21(«)-t-j  -«-Z2-t  ,,(«.) 

z_- and       z_(«).= 


velocity 


Z2t   I2(<l))-(-j      tt)t  22((D) 
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In  order  to  determine  a  numerical  solution  for  the  analysis,  physical  characteristics 
of  the  system  must  be  input.  These  values  used  here  and  throughout  the  thesis  are  based 
upon  a  real  world  oceanographic  mooring  described  by  Vandiver,  discussed  in  chapter  4. 

The  following  are  the  physical  characteristics  of  this  system: 


a    =  .0254  ra 


b   =.0254m 


p  .    =1.144^ 
m 


d=.1524m  c:=.  1778  m 


f=.1524m 


e  =.  1778m 


h    =c-(-d 


i:=e+ 


M  j  :=4.52  1c 


M2:=4.52k 


p2=1.144-lg  Z,:=33.8231-^         Z - :  =33.1231  JS        K=10^^ 

m  sec  sec  m 


M  ,(2d)" 


12 


R    =0- 


M2(2f) 

12 
newton 


H    =100Oncwton 


Pi 


P2 


A  curve  of  the  Power  and  Transmission  Coefficients  provides  an  indication  where  to 
search  for  the  Resonance  Frequencies  of  the  system 


100 


150 


200 


250 


(radians/ sec) 

Figure  A-2:  Graph  of  Reflection  and  Transmission  Coefficients  Using  Input 

Parameters  from  Table  4-1  (R=0) 

Based  upon  this,  it  appears  as  if  the  best  absorption  would  occur  at  the  frequencies 
at  the  location  of  zeros  of  the  Power  Reflection  Coefficient.  Using  the  secant  method,  we 
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try  to  obtain  these  values  at  the  location  of  the  curves  where  the  Power  Reflection 
Coefficient  appears  to  go  to  zero. 

The  reflection  coefficient  can  be  reduced  to: 


24222222222  32 

l-t2,((i))   +(i)    Z  ,    Z-,   t  )2(u)   +  (0    Z  ,    t22<<i>)   -  0)    Z2   t  j  ,(o>)  +  2ruZ  ,  t  ^^cu»t  ,,((o  )  +  2r(D    Z  ,  Z2   t  t  ,(<i>)  t  )2(<i>) 

RclCi))     = - — ! 

2  2  42222222  22  2 

t2](<i>)   -  2  l-,jt(D  i  (i)    Z  ,  Z2t  j 2< (!))  +  <■>    Z  j    Z-,    t  ]t((i))   +u    Z  ,    t  22((i))   +2  (i)    Z  j  t  22d!))Z-,t  ,|((i))+(i)    Z  ,    t  ||((i)| 


And  to  determine  the  zeros,  we  solve  for  the  values  of  w  such  that  the  numerator 
goes  to  zero  using  the  secant  method. 

The  first  data  point  is  close  to  40  rad/sec  for  the  initial  guess  <■>  =42.564  — 


The  expression  for  the  solution  is  provided: 


"iRc 


root 


2422  222  222  2\ 

t  2i(u  )    ■+■»    Z  j    Z2    t  12(<i))    -t-o  -Z  j    t  22<w)    -  «>    Z2    I  n(«>)  j  -l-2i(i)-Z  j  t  22(0))'t  2l(C0)  -  •lu 

+  2iw   -Z  j  Z2    t  n(a))t  ,2(u) 

The  resulting  resonance  frequency  is:  »  =42.564  — 

sec 

This  convergence  represents  an  iteration  between  estimate  and  calculation  until  the 
solution  converged  to  the  third  decimal  point. 

The  second  data  point  is  close  to  90  rad/sec  for  the  initial  guess        &>  =94.13  — 

sec 

The  expression  for  the  solution  is  provided: 


Q2Rc 


root 


2422  222  222  2\ 

t  21(o> )    -^-0)    Z  y    Z  2    t  12(0))    -Hw    Z  j    t  22(w  )    -0)  -Z2    t  n((i))    I  -l-2i-(i>-Z  j  t  22( «  )t  21^w  ^  •■■■.<* 

3  2 

+  2iw    Z  j  Z2    t  n(«)t  ,2((o) 


The  resulting  resonance  frequency  is:  <o  =94.13 


rad 


This  convergence  represents  an  iteration  between  estimate  and  calculation  until  the 
solution  converged  to  the  third  decimal  point. 
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A  second  determination  of  the  natural  frequencies  relies  on  the  determination  of 
the  value  of  w  when  the  input  impedance  matches  that  of  the  input  string.  Because  the 
string  has  a  purely  real  impedance,  the  value  of  co  such  that  the  imaginary  part  of  the  input 
impedance  is  zero  will  provide  this  frequency. 

The  imaginary  part  of  the  input  impedance  is  described  by: 

3,    2 
-co  >Zo  't  n((i>)'t  12(0)  — I  2j(°  )'ffi't  ">2^<l)  ^ 


^(imaginary)  =  ■4-z22-t1rf->2+.2-ta<.>2 

Again,  the  basis  for  the  determinations  of  the  value  of  co  such  that  the  expression  is 
zero  is  the  location  of  the  minimum  values  of  the  Reflection  coefficient  observed  from  the 
graph,  and  using  the  secant  method,  we  can  determine  the  resonance  frequencies. 

The  first  data  point  is  close  to  40  rad/sec:       a.  =42.564 


rad 
sec 


The  expression  for  the  solution  is  provided: 


fl  1  Zm    ~rootiv"(l)    z2    l  H*®)"1  12(<")-t21((i>)-(0-t  22<(0)-(1" 


The  resulting  resonance  frequency  is:  s>  =42.564 


rad 


This  convergence  represents  an  iteration  between  estimate  and  calculation  until  the 
solution  converged  to  the  third  decimal  point. 


The  second  value  appears  to  be  around  90  rad/sec 


The  expression  for  the  solution  is  provided: 


u    =  94.13 


rad 


(     3       2  \ 

A  2  Zm    _root'_u    z  2    '  n* u  H  12(M  )  ~  l  2l(°)  *'0)'t  22(u  ),tB  ; 


The  resulting  resonance  frequency  is:  <u  =94.13  — 


rad 
sec 
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This  convergence  represents  an  iteration  between  guess  and  calculation  until  the 
solution  converged  to  the  third  decimal  point.  For  the  impedance  method,  the  result  was 
not  as  sensitive  to  the  initial  estimate  of  the  value  of  the  frequency,  but  the  Reflection 
Coefficient  and  the  Impedance  Method  provided  the  same  results. 

A  sketch  of  the  two  mode  shapes  that  correspond  to  the  location  of  the  reflection 
zeros  are  provided  in  Figure  A-3.  The  sketch  represents  an  estimate  of  how  the  two  mode 
shapes  can  be  visualized,  although  a  solution  to  the  total  system  motions  was  not 
performed  in  order  to  provide  this  sketch.  Position  A  and  position  B  are  the  two  extremes 
of  the  mode  shapes.  The  slight  relative  rotation  between  the  two  masses  in  mode  1 
represents  the  amount  of  absorption  that  occurs  in  this  mode  shape.  The  extreme  relative 
rotation  between  the  two  masses  in  mode  2  represent  the  maximum  amount  of  absorption 
that  can  take  place  in  the  system  at  the  second  resonance  frequency. 


Cable 


Cable 


Mode  1 :  Slight  Absorption 


Cable 


Cable 


Mode  2:  Maximum  Absorption 

Figure  A-3:   Mode  Shape  Sketch  of  System 
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Appendix  B 

Determination  of  the  Optimum  Damper  Values 

Based  upon  the  formulations  of  the  transfer  matrix  equations  from  Li,  the 
following  analysis  describes  the  approach  to  determining  the  optimal  value  of  the  damper 
operating  at  the  second  resonance  frequency.  A  computer  program,  Mathcad  5.0®,  was 
used  in  the  calculation  process.  The  first  difference  between  this  formulation  and 
Appendix  A  is  the  fact  that  the  relationships  are  dependent  upon  R,  not  co. 


a  1<R>    =T  i  minus  d  +  Hc+Ka2-(w)    I  ,  -hj    wRb2 


2  2  2 

a  2(R)   ~T  jf-l-He  +  Ka  -(&))    I2+-j    ffl-R-b 


a3(R)  :=K-»2+j    wRb2 

Intermediate  expressions  are  described: 


T  i(R)=- 


(M  jl-t-Mvfja  j(R)-M  j-c-(a3(R)W<fl    M  ,M  2fc2 
M2-e-(a  ,(R))  -t-(M  j-d-1-M2-h)(a3(R))-f-(0    M  ,M  2<?c-d 


,X-«,..a^^ 


T  2<R>:=- 


M2-e-(a  j(R))-(-(M  |-dfM2'h)-a3(R)  H-(i)    M  ,M  2c<rd 


M  ,-dlc+l 


(a  1(R))-l+(a3(R)  -h 


T  3<R> 


j  M2-c-(a  j(R))  +  (M  j-d-(-M2-hWa  3(R))  H-o>    M  j-M  2ced 


T  4(R) 


(M  j-l  +  Mvfj-(a3(R)J  4-M  [■ca2(R)-l-(ii    M  j-M  2e-c-f 


OjIRVh  +  OjtR)-! 


M  j-d-^+M^-h 
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The  intermediate  expressions  continue: 


Y  s(R>:=- 


(a2(R)h  +  a3(R)l) 


+  M  ye-fh 


o^RVh  +  a^RVl)      Lm        2      m    ^x 


M  j-d-l    -hM2rhl 


T6<R)    = 


2  2 

-a^RJ^M  jd-(-M2h)-M2e-a3(R)-(-iB    M  ,M  2<s    d 


(tt2<R)h  +  a3(R)l) 


M  .•d-l2fMyf2'h 


Now  knowing  these,  the  transfer  matrix  can  be  determined: 


T  5(R)+Y  2<r)Y  6<R> 
1-7  3(R)T  6(R) 
7  4(R)  +  Y  i(R)T  6(R> 


1  12<R>    ~T  2<R>+T  3<R) 
t2l(R) 


1-7  3(R)Y  6(R) 

Now  the  reflection  coefficient  is  determined 

_4(t21(R))"t"(B    Zl'Z2t  12<R)   "t-j  •»(Zr'22(R)-Z2'lll(R)) 


t  n(R)  -T  !(R)+7  3(RV 
l22<R> 


T  4(R)+7  ,(R)Y  6(R) 


1-7  3(R)7  6(R) 
7  5(R)  +  7   2<R)Y  6(R) 


1_T  3(R)Y  6(R) 


RC(R)=- 


t21(R)-o)     Z  ,Z2t  12(R)+j    u  (Z  ,t22(R)+Z2t  n(R)) 

The  transmission  coefficient  is  determined: 


TC(R):=- 


j  -2-b-Z 


t21(R)-a     Z  j  Z2t  12(R)-(-j    »(Z  jt22(R)+Z2t  n(R)) 


The  impedance  is: 


zm(R) 


t2,(R)+j    w  Z2t  n(R) 

2 
-to   -Z2t  12(R)-l-j    wt22(R) 


The  Power  Reflection,  Transmission  and  Absorption  Ratios  are  calculated 


P(R):=RC(R)-RC(R) 


7(R)  •=—  TC(R)TC(R) 
Z, 


a(R):=l-p(R)-y(R) 
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Now  the  determination  of  the  optimum  value  of  the  damper  value  to  obtain  the  most 
absorption  can  be  done  by  taking  the  derivative  of  the  Power  Absorption  Ratio  with 
respect  to  R. 


knowing 


-RCR 


or 


5a   _    ^p 6y_ 

5R  5r     6r 


_Z2        - 
T- — T    -T 


Because 


a      =     1-P-t 


leads  to: 


R  CR  J 


5a 

5R  5R 


5R 


5a     _        «Rc—  5RC       Z2    5TC_  5T  J 

—     ~ Rc+Rc T.tT. 

5R  \5R  5R  /      Zl    \5R  5R  / 


This  derivative  requires  the  calculation  of  the  derivative  of  the  transfer  matrix  from  the 
constituent  components  of  the  system.  The  intermediate  equations  are  provided. 


1  l_prime_«(R>=[Mr<I+c)+M2f]'  M2-e-(a  j(R))  +  (M  ,d  4-M  2-h)-  (a  3(R))  +  o>2-M  rM  2c-e-d 
1  l_prime_b(R>  ;=[M  i'd  +  M2(e+h)]|-(M  jl-HM^fj-a  ,(  R)  -  M  ,-c-  (a  3(R))  +  a>2-M  ,-M  2f  c 
J  ■••^■(t    l_prun._,(R>+lf   l_prune_b<R>) 


*   l_primc(R)    : 


M2-e-ia  j(R)j  +  (M  j-d-f-M  2'hWa  3(R)J  +w    M  j  M  2-c-e-d 


T2_prin,e_.(R>:=-   ^ '•[  M  2^  (a  ,(  R)j  +  (m  ,-d  +  M  2h)- (a  3(  R)) +o)2  M  ,-M  2-c-cd  | 
\  o)     / 


M_pnme_b(R):=[Mld  +  M2-(e+h)] 


-M2fh  -M  ,-l-d   + 


[(a,(R)Vl+(a3(R))-h]' 


2_primc 


(R) 


J    fflb     (t  2_pnme_a(R^+l'  2_pnmeJ><R>) 


M2-c-(a  ,(R)j  +  (M  ,-d-(-M2-hJ^a3(R)j+(i)    M  jM2ce-d| 
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T3_primc_.(R);  =  (— ]-[M2-c(o1(R))-(-(M1-d  +  M2-h)-(o3(R))-t-»    M  ,-M  2c-cd 


T  3_prunc_b<R>    =[M  ,d +  M  2(e  +  h)  I 


(a  1(R))-l+(a3(R))-h] 


1  3_primc(R) 


"J  ^-{l  3_pnmc.a(R>^Tr  3_pnn,eJ,(R>) 


l[-M2-e-(a  i(R))-  (M  1-d  +  M2-h)-(a3(R))]-(o    M  ,M  2ce-d 


1  4_prime  a(R)    : 


[M  j(l-t-c)  +  M2fl 


(a2(R)h-f-a3(R)l) 


-  (m  (d-^  +  M^^h 


4_primc_b 


(R) 


2   j 


(M  j-H-M2-^-(a3(R))+M  j-c-a2(R)+(i)    M  j-M^e-c-f 


4_prune 


(R) 


"J    wb     (T  4_pnmc_a(R)+T  4_pnmcJ><R>) 


(a2(R)-h-i-a3(R)-l) 


VI  .d^  +  Myf^h 


.2 
j    Mb  • 


5_primc 


(R) 


h+,  r(a2(R)-h+a3(R)-l) 


-  (m  j-d-^  +  Mj-^-h 


(a2(R)-h+a3(R)IJ 


-t-M2e-fh 


]J 


(a2(R)h-ha3(R)l) 


IM  i-d-^  +  Mv^-h 


T  6_prime  a(R) 


-[M  j-d-(-M2-(ei-h)l- 


(o2<R)-h  +  o,(R)-l)       /  ,  , 

- --    M  ,-d-l    fMvPh 


iU  I     -p  ivi  2 


6_prime_b 


(R>:  = 


h  +  l 
2 

0)      /J 


2  1 

■a  2(R)-(M  jd  +  M  2-hj  -  M  -,-c-a  3(R)+-w    M  ,M  2e    d 


1  6_pnme(R 


j  '"'•^'(t  6_pnmc_a(R)+T  6_pnme_b(  R>) 


(o2(R)-h-t-a3(R)-l)       /  ,  y 

(m  j-d-i  +M2r-h 
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The  elements  of  the  derivative  of  the  transfer  matrix  are  then  calculated: 


^Ijmmc.,^)^  4_pnmc(R>  +  T   ijaimeW'T  6<R>+Y   1<  R>T  6_pnmc(  R>  +  T  3_prime<R>T  4<R>  T  6<  R> 

t21_pnme.b(R)  :=-T  3<R)T  4_prime<RH  6<R>  +  T  3<RH  4<RW  6jrime(R>    ' 
+  (T   1<RH  SjW^-T   l_primc(R)T  3<R>H*  6<R>)2 

.'21jnBi«  i(R)+t21j)nnieb(R) 


(R) 


21_prime 

(l-T  3(R)'T  6(R>)* 

t22_primc_.(R):=r  5jnnK(R)  +  T  2_pnme(RH  6(R>+T  2<RH  6_prime<R> +T  3_pnme(RH  5<RH  6(R) 

t22_pnme_b(R)=-l'  3<R>*   5_prime<R>-T  6<R>  +  T  3<RH   5<RH  6_pnmc(R>" 
+  (T  2<R>?  3_pnme(R>-T  2JJrin,e(R)T  3<R>)'(r  6<R>)2 

_ l  22_priroc_i(  R)  "Ht  22j>rime  b(  R) 


t22_prime(R 


l-T  3(R)-T  6(R)J 


*  ll_prime<R>  :=^   lj»une(R)  +  Y  Sj^nc^H  21<R>+T  3<R>'  21_pn«ne(R> 
l12_primc(R);^  2_prim*(R^T  3_pnmc(R)',22(R)^T  3<R>  <  22_prunc(R> 

The  calculation  of  the  derivatives  of  the  Reflection  and  Transmission  Coefficients  with 
respect  to  R  are  provided: 

dRc_dR>R)  :=<D2-Z2-(t  ^R)-t21jw^(RH^ 


dRc  dR  WR)  :=j  • 


3,2/ 


dRc_dR(R) 


"    Z2     (tnjmmC(RVtl2(R)-hl(R)h2J>rinM(R))-ffl(t21^me(R)l22<R)-t2l(R)t22_prin,e(R)) 
2Z  ,(dRc  dR  a(R)-f-dRc_dR_b(R))  SRc 


t21(R)-w    Z  j-Z2t  J2(R)I  -fj    (.)(t22(R)Z  j  -t- 1  n(R)Z2] 


6R 


dTc_dR(R)    =2Z  ,■ 


M   •(t22pnme(R)Zl  +  tll  pr»nc(R>Z2)+J  >    Z  1  Z  2  '  12  pn,nc<R>  ~  »  '  21  primc<R> 
[(t21(R)-ffl   -Z  ,Z2t  12(R))+j    (u(t22(R)Z  j  -Ht  U(R)Z2) 


_       STc 
5R 


The  operation  of  the  complex  conjugate  and  the  partial  derivatives  derivative  are 
interchangeable,  which  leads  to: 


5a  _    /. 


2  /, 


—  =  -ldRc  dR(R)RJR)+R  „(R)dRc  dR(R) dTc  dR(R)T  r(R) -(-T  „(R)dTc  dR(R) 

6R        V       "  '      Z,    ^       -  c 
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Again,  in  order  to  complete  a  numerical  analysis,  a  set  of  physical  characteristics  of 
the  system  must  be  input.   In  the  case  of  this  example,  the  previous  values  from  Table  4-1 
are  used  again  with  a  concentration  at  the  second  resonance  frequency  of  the  baseline 
system.  Evaluated  at  the  higher  of  the  two  resonance  frequencies: 


a   =.0254  m 


b    =.0254m 


p  {    =  1.144-^ 
in 


d=.1524m 


f=. 1524m 


p  2    =  1.144- 


kg 


c   =.  1778  m 


e:=. 1778  m 


h    =c  +  d 


I    =e  + 


M  j    =4.52k 


M2=4.52lc 


ke  lea  newton 

Z.    =33.8231  —        Z  2   =33.8231—         K=10 

sec  sec  m 


2    " 


M  j  (2  d)^ 

12 

M2-(2f)2 

12 


nd 

0)    =94.13 


H   =100On<rwU)o 


Pi 


t  ;  :=- 


p2 


In  order  to  obtain  an  initial  estimate  of  the  value  of  the  damper  for  maximum 
absorption,  a  graph  of  the  Power  Adsorption  Ratio  is  provided.  This  will  provide  the 
range  for  the  initial  estimate 


400  600 

R(kg/»ec) 


1000 


Figure  B  -  1:   Plot  of  Absorption  Power  Ratio  For  Varying  Values  of  R  Using  The 

Input  Characteristics  of  Table  4-1 
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The  expression  to  determine  the  optimum  value  for  the  damper  is  provided  by: 
The  initial  estimate  is  around  500  newton/(m/sec)  from  the  graph  above 


Bart  R 


R    =464.343 

m 

sec 

root 

-  ^dRc_dR(R)R  C(R)  +  R  c(  R)dRc_dR(  R)j  ... 

.R 

+  - 

Z  2  / 

J 

—  (dTc_dR(R)T  C(R)  +  T  c(R)dTc_dR(R)j 

newtaa 

Bert  R  =464.343  *— — — 

This  represents  an  iteration  until  the  result  was  within  three  decimal  places  of  the  estimate. 
An  iteration  was  performed  until  the  solution  was  within  three  decimal  places  of 

the  initial  estimate.  The  value  of  the  optimum  damper  for  this  case  is:  464.343  kg/sec. 
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